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Abstract. Extensions of the classical orthogonal polynomials which satisfy 
differential equations have a long history. Recently, a new soliton method 
was used to construct commutative algebras of differential operators having 
orthogonal polynomials as eigenfunctions. The goal of the present paper is to 
extend these techniques to the Askey- Wilson level and to describe explicitly a 
large collection of families of orthogonal polynomials which are eigenfunctions 
of commutative algebras of g-difference operators. 



1. Introduction 

Orthogonal polynomials which are eigenfunctions of a differential operator have a 
long history. In 1929, Bochner ^ proved that, up to an affine change of variables, 
the classical orthogonal polynomials are the only ones that satisfy a differential 
equation of the form 

B ^X,-^^ Pn{x) ^ \nPn{x), n G Nq, (1.1) 

where B{x,j^) — J2j=o ^ second-order differential operator with coeffi- 

cients bj{x) independent of the degree index n. In 1938, Krall [17] posed the general 
problem to construct and classify all families of orthogonal polynomials which are 
eigenfunctions of a differential operator B (a;, ^) of arbitrary order, which is inde- 
pendent of the degree index n. He showed that the order of the operator must be 
even and solved it completely for operators of order four [TH]. During the next 60 
years, many different examples of orthogonal polynomials satisfying higher-order 
differential equations were constructed, see for instance [Mj [TBI [20] and the refer- 
ences therein. 

The pioneering work of Duistermaat and Griinbaum brought new bispectral 
techniques in this old problem. Roughly speaking, the bispectral problem concerns 
the construction and the characterization of functions depending on two variables 
which satisfy simultaneously spectral equations in both variables. For the classical 
or Krall polynomials, the bispectral property can be explained as follows. Recall 
that for any family of orthogonal polynomials, the operator multiplication by x can 
be represented by a three-term recurrence relation: 

XPn{x) = AnPn+lix) + BnPn{x) + C„p„_i(x), (1.2) 

where the coefficients An, Bn, C„ are independent of x. Looking now at equation 
(|l.ip we see that the polynomials are eigenfunctions of a differential operator acting 
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on the variable x, while equation ()1.2[) shows that they are also eigenfunctions of a 
second-order difference operator acting on the degree index n. 

Although Duistermaat and Griinbaum considered differential operators in both 
equations, the connection between Krall's problem and the bispectral problem, 
suggested that the tools used in [5] can be translated within the context of orthog- 
onal polynomials. And indeed, the general bispectral techniques developed in [3] 
were used to construct far reaching extensions of the Laguerre polynomials [7] and 
Jacobi polynomials which are eigenfunctions of higher-order differential opera- 
tors. Things were taken even further in [5] where extensions of the Askey- Wilson 
polynomials which are eigenfunctions of g-difference operators were constructed, 
containing as g — > 1 the previous results. The core of the all these papers is the 
so-called Darboux transformation, one of the basic tools in the theory of solitons. 

It is interesting that yet another soliton method was efficiently used to con- 
struct bispectral operators. It appeared first within the context of general Dar- 
boux transformations for differential operators [22], and a few years later it was 
more systematically used to prove the bispectrality of specific rational solutions of 
the Kadomtsev-Petviashvili hierarchy [19] and Toda lattice equations [10]. It was 
adapted very recently in [TlJ [12] to establish the bispectrality of generalized Jacobi 
and Laguerre polynomials. One of the advantages of this approach is that it pro- 
vides a more detailed information about the commutative algebras of differential 
operators which was crucial for multivariate extensions |11| . Moreover, it naturally 
leads to the construction of lower order differential operators, thus confirming con- 
jectures concerning the characterization of the commutative algebra of all possible 
differential operators diagonalized by a specific family of orthogonal polynomials 

The goal of the present paper is to extend the techniques from [TTl [12] to the 
Askey- Wilson level. In particular, this provides a new approach and different proofs 
of the results in [9] . It also explains from a general point of view the existence of 
some lower-order g-difference operators found there in a rather roundabout way 
for specific examples. Moreover, we have organized the presentation in a slightly 
more general setting, in order to include also polynomials satisfying higher-order 
recurrence relations (thus, in particular, the constructions here include polynomials 
orthogonal with respect to Sobolev-type inner products). Since the Askey- Wilson 
polynomials are at the very top of the g- Askey scheme, the statements in the present 
paper lead to a variety of interesting special or limiting subfamilies of orthogonal 
polynomials which are eigenfunctions of q-difference, difference or differential oper- 
ators. Producing an explicit table with all such examples similar to the reductions 
of the Askey- Wilson polynomials [15] is a challenging task, far beyond the scope 
of this paper. On contrary, we have tried to make the presentation compact and 
self-contained, emphasizing the main new ideas and constructions. 

The paper is organized as follows. In the next section, we fix the notation 
and collect the main properties of the Askey- Wilson polynomials needed in the 
paper. In Section [3] we construct the extensions {Pnix)}^^Q of the Askey- Wilson 
polynomials and formulate the main result - their bispectrality. In Section |4] we 
establish the recurrence relations (the spectral equations in the degree index n). 
The techniques here are more or less known, so we outline the proof by stressing 
the specific points and how they relate to the existing literature. Section [5] is 
the heart of the paper where we construct the commutative algebra of (/-difference 
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operators in z, with x — ^(^ + which is also diagonahzed by the the polynomials 
{Pn{x)}^^Q- Finally, in Section[B]we treat in a detail extensions of the Askey- Wilson 
polynomials orthogonal with respect to a measure on the real line. In particular, 
we obtain a different derivation and a refinement of the results in [9J. We also 
formulate the main statements in terms of orthogonality. 



2. Notations 



Throughout the paper we assume that q is a real number in the open interval 
(0, 1) and we consider the corresponding g-shifted factorials 



{a;q)n = ]^(l-ag'), {a;q)oo = JJll-a?') and (ai, 02, . . . , ; g), 

1=0 1=0 

and 41/13 basic hypergeometric series 



493 



ai, 02, 03, 04 
61,^2,63 



q,z 



E°° (ai, 02, as, 04; (j); ; 
- ' ' —z . 



1=0 



{bi,b2,b3,q; q)i 



The Askey- Wilson polynomials Pn{x) = Pn{x; a, &, c, d), n ~ 0,1,2.. 
in [2] depend on four free parameters a, 6, c, d and can be defined by 

{ab, ac, ad; q)r 



p„(x;a, 6,c, d) 



-493 



q ",abcdq" ^,az,az ^ 
ah, ac, ad 



;q,q 



Y[{ai;q)n, 

1=1 

(2.1) 
introduced 

(2.2) 



where x — i(z + i). Sears' transformation formula [SI page 49, formula (2.10.4)] 
shows that the polynomials Pn{x;a,b,c,d) are symmetric in the four parameters 
a,b,c,d. We denote by £"„ and E~^, respectively, the customary forward and 
backward shift operators, acting on a function /„ — f(n) by 

If we define the second-order difference operator 



Cn — AnEn + -BrJd + CnE^ ^, 

with coefRcients 

1 - ahcdq^-'^ 



Br, 



(1 - a6cdg2«-i)(i _ abcdq^'^) 

(7"-i[(l + abcdq^''-^){sq -t- s'abcd) - g"-i(l + q)abcd{s + s'q)] 
(1 - a6cd(72"-2)(l _ abcdq^'') 



(2.3) 



(2.4a) 



(2.4b) 



Cn = (1 - a6g"-i)(l - acg"-i)(l - adg"-i)(l - 6cg"-i)(l - 6rfg"-i)(l - cdg"-!) 



(1 



(1 - abcdq'^"--'^){l - abcdq^"^-^) 
-1 I I . 



(2.4c) 



where s = a + b + c + d, s' = a +6 



then we have 



^nPnix) = 2xpn{x). 



(2.5) 

Note that the action of £„ is well defined on functions on No (since Cq = 0). It 
will be convenient later to use a natural bi-infinite extension of C„ which acts 
on functions defined on Z. We can do this by the formal change of variables 
n— ^7 = n-|-ein the coefficients of Cn ■ We shall consider e for which the coefficients 
A^, i3-y, C-y are well-defined (i.e. the denominators do not vanish) and A-^ ^ 0, 
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C-y ^ for n e Z. We shall denote the bi- infinite extension of £„ by C-y. We use 
the same convention for other difference operators acting on n, i.e. if 

finitely many j 

is a difference operator acting on functions defined on No, we denote by 

finitely many j 

its bi-infinite extension acting on functions defined on Z, where 7 = n + e and 

Ej = En. 

Let us denote by and Dj^, respectively, the forward and backward g-shift 
operators, acting on a function h{z) by 

Dzh{z) = h{qz) and D^^h{z) ~ h{z/q), 

and let Bz be the Askey- Wilson second-order g-difference operator 

= A{z)D, ~ [A{z) + A{l/z)] Id + A{l/z)D-\ (2.6) 



^ _ (J- - - bz){l - cz)il - dz) 

(1 — z^){l — qz^) 



BzPn{x) = XnPn{x), (2.8) 

= (g-"-l)(l-a6cdg"~i). (2.9) 



where 

Then 
where 
We also set 

^ q''iabcd/q;q)n ^ q''{abcd/q, hcq^ M(f\ cdq"^, g)oo 
" a^'-{bc,bd,cd,q;q)n a'^{abcdq^~^ ,bc,bd,cd,q;q)Qo 

The first formula will be sufficient for most applications; the second representation 
will be used when we want to consider a meromorphic extension of in n (and, 
in particular, when we want to make the change of variable n — > 7). Note that 
is symmetric in b,c,d and we denote by the analogous expressions, 

where the parameters (a, b, c, d) are replaced by (6, c, d, a), (c, d, a, 6) and (d, a, 6, c), 
respectively. 

Using the 4(^)3 representation of Pn{x) given in (12. 2p . one can show that the 
polynomials ^^p„(x) satisfy the following q-difference-difference equation 

B'^M>^{x)-Cn-iPn-i{x)) - (q-"+i-a6cdg«-i)(e>„(x)+dp„_i(x)), (2.11) 

where 

"^^'^ -A''{z)-A''{l/z) 



Bl = A''{z)Dz + 
and 

A^z) 



q 

{q + az){l - bz){l - cz){l - dz) 



Id + A-il/z)D;\ (2.12) 



(1 -z2)(l -g^2) 

We denote by B^, B^, Bf the operators analogous to B'^, where the parameters 
{a,b,c,d) are replaced by {b,c,d,a), {c,d,a,b) and {d,a,b,c), respectively. Thus 
equation (|2.1ip holds for any permutation of the parameters (a, b, c, d). 
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We use Wr„ to denote the discrete Wronskian (or Casorati determinant) of func- 
tions hn \ hn \ . . . , h!"n^ of a discrete variable n, i.e. 

Wv^{h(^\h[^\ /iW) = det(/iiVi)i<,,,<fe. 
3. Statement of the main result 

Consider now arbitrary polynomials (j)^^"^ (A„) € C[A„], for j — 1,2, . . . ,k, and for 
each j pick Sj G {a, b, c, d}. We set 

= (3.1) 

We shall assume that ipn^ are independent (as functions of n). We define a new 
family of polynomials in x by the Wronskian formula 

Pn{x) = Xn Wr„(V'(i),7A(2)^ . . .,^i^\p^{x)), (3.2a) 

where x„ is an appropriate normalizing factor (depending only on n). The main 
result is that the polynomials Pnix) are simultaneously diagonalized by two com- 
mutative algebras. The first algebra consists of difference operators acting on the 
variable n (with coefficients independent of z), while the second algebra consists 
of g-difference operators acting on the variable z (with coefficients independent of 
n), where x = \{z -\- \), thus extending the bispectral equations (|2.5I) and (|2.8p 
satisfied by the Askey- Wilson polynomials. The factor Xn is not essential for the 
bispectrality, but we shall fix it below in order to obtain well-defined polynomials 
for all n G No and to formulate the precise statement. 

First, note that if we multiply the j-th row in the determinant on the right-hand 
side of equation (|3.2a[) by Cn-j+i? for j = 1, . . . , fc + 1 we can rewrite the formula 
for pn {x) as follows 

k+l , 

Pn{x) = Xn n 7S Wr„(CV'i'),e>f (3.2b) 

Let T„ denote the principal fc x fc sub-determinant 
fc+i , 

= Xn n ^ ' ^ ' • ■ • ' ) • (3-3) 

If 5j — a, then from (|3.ip we see that in the j-th column in the above formula we 
have Laurent polynomials of g". If j = 6 wc can pull in front of the determinant 
the factor and the remaining entries will be rational functions of g". 

It is easy to see that if we pick 

= ^;iT) iacq--''^\adq-^+^;q),.„ (3.4) 

then a multiplication by 7/^'^ cancels all denominators in the j-th column and the 
entries become Laurent polynomials of g". We define ri!f^^ and yy^'-' by replacing in 
p.4p the parameters {a,b,c,d) with {a,c,b,d) and {a,d,c,b), respectively (i.e. we 
keep a fixed, and we exchange the roles of b and c, or b and d). Thus, if we define 

k+l 

Xn^Yl^tj+i n (3.5) 

j=l l<3<k 
Sj 
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where we set Sk+i — a, then t„ becomes a Laurent polynomials of g". 
Throughout the paper we shall assume that 

r„ 7^ for n = —1, 0,1,..., (3.6a) 

which is true for generic polynomials 0'^-')(A„), and 

^"^"-''^^ ^ for n = 0, 1, . . . , (3.6b) 

Xn-l 

which is true for generic parameters (a, &, c, d). 

Next, we define two algebras 21^ and 2l„ of Laurent polynomials in z and g", 
respectively. Let 

W^ = span{V;W,Vf (3.7) 

be the space spanned by the functions V'V'' ^^"^ denote by 2tz the subalgebra 
of C[z + 1/z] consisting of all polynomials f{z + 1/z), for which the space W is 
/(£-y) -invariant, i.e. 

21, = {f{z + 1/z) e C[z + 1/z] : f{C^)W C W}. (3.8) 

The fact that 21, contains a polynomial of every sufhciently large degree follows 
from Proposition 14.11 Interesting examples of the above construction when 21, = 
C[z-|-l/z] and the polynomials {pn{x)}'^^Q are orthogonal with respect to a measure 
on R are discussed in Section [HI 

Next we define the algebra 2l„ which consists of all polynomials ft.(A„_j./2) G 
C[Xn-k/2] such that h{X,-,_k/2) - H\i-k/2-i) is divisible by t„_i in C[g",g~"], i.e. 

2l„ = {mA„-./.) e C[A„_,/,] : e Cr,g-"]} . (3.9) 

The fact that 2lri contains a polynomial of every sufficiently large degree follows 
from Proposition 15. II 

The main result of the paper is the following theorem. 

Theorem 3.1. For every f{z + 1/z) G 21, and every h(Xn-k/2) G 2l„ there exist 
a difference operator acting on n, which is independent of z, and a q-difference 
operators'^ acting on z, which is independent of n, such that 

^iPn{x) = f{z + l/z)p„{x), (3.10a) 

B'^Pnix) = h{X,_k/2)Pn{x). (3.10b) 

Thus, Tin = {C-li^ : / G 21,} is a commutative algebra of difference operators in n 
isomorphic to 21,, S), — {B'^ : h e 2t„} is a commutative algebra of q-difference op- 
erators in z isomorphic to 2l„, and both algebras are diagonalized by the polynomials 
Pn{x). 

We establish equations p.lOa|) and (j3.10bl) in Sections H] and [SI respectively. 

4. The difference equations in n 

Consider the bi-infinite extension of the Askey- Wilson recurrence operator 
defined by equations (|2.3p - (|2.4p and for each £ e {a, 6,c, d} we denote by the 
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difference operator obtained by conjugating with the operator multiphcation by 
i.e. 

>C^=C^/:1. (4.1) 
Exphcitly, when f = a we can write as follows 

£; = ^^S7 + B;id + c^i;-\ (4.2) 

where 

^ a(l - 6cg^)(l - „ cdg^)(l - g^+i) 

g(l-a5cdg27-i)(l_a6cdg27) ' ""^ ^ 

= - + - - - (4.3b) 

^ g(l - ahcdq^-^)(\ - abg^-i)(l - acq^-^)(\ - adq^-^) 
^ a(l - abcdq^^-^){l - ahcdq^^-^) ' ^ ' 

From these formulas it is not hard to see that for every j e No we have 

{X-yY = {—pfi ^ ) ('^7)"' + ^ polynomial in of degree at most (j — 1), 

where A-y is the eigenvalue of the operator Bz defined in (|2.9p . The above equation 
shows that every element in C[A7] is annihilated by an appropriated polynomial of 
£°. Combining this with p.ip . (13. 7p and (14. ip we obtain the following statement. 

Proposition 4.1. There exists a nonzero polynomial f such that f(C^)W = {0}. 
In particular, the algebra 21^ defined by p.Sp contains a polynomial of every suffi- 
ciently large degree. 

We can easily write an explicit eigenvector for the operator £° in the space C[Xj] 
for the eigenvalue a/q'^'^ + q^~^-^ /a. 

Proposition 4.2. For j E No define 

'q-^ , q^+^/a^, q'^+\l /abcdq''-^ 



g]i\) = (q /ab,q^/ac,q^/ad;q)j i(l33 



qyab,q^/ac,qyad ''^''^ 



(4.4) 

Then 

^7ff"(^7)=(^ + ^)ff"(A7). (4.5) 

Proof. Using the definition (j2.1l) of the 4(^3 basic hypergeometric series, one can 
check directly that (|4.5p holds. □ 

Note that if a ^ {ig*/^ • g ^ Nq}, then g^iX^) is a polynomial of degree j in A^. 
In this case, the eigenvalues a/q^^^+q^^^ /a are different for j e Nq, the operator 
can be diagonalized in C[A'y], and the polynomials given in ()4.4p provide an explicit 
eigenbasis. Clearly, we obtain similar statements for the operators >C5j,, C^, £ij and 
we denote by ^^(A^), ^^(A^), g^(A-y), respectively the corresponding polynomials in 
formula (14.41). 



Remark 4.3. We note that the polynomial gj{Xj) can also be obtained from 
[13] . Indeed, appropriately normalized, the polynomial gjiXj) in equation ()4.4p 
above corresponds to the eigenfunction Sf{aq^^^^) of the bi- infinite Askey- Wilson 
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recurrence operator given in formula (1-13) on page 203 in [13]. To see this, we 
need to multiply s-y{aq~^^^) by an appropriate factor depending on 7 (because 
the recurrence operator there differs from £°). Then we apply first the Watson's 
transformation formula 6, formula III. 18] to rewrite the 8(^)7 series as series, 
and then we need to use Sears' formula [6j formula III. 16] to transform the 4^3 
series. 

Consider now the difference operator Qj defined by 

Q7/7 = X7 Wr^ (V^') , . . . , V-^'^ , A). (4.6) 
Note that W = ker(Q^) and therefore, for / € 2lz we have 

/(£^)ker(Q^) Cker(Q^), 

which implies that 

ker(Q^) c kei{Q^f {£-,)). 

This means that there exists a difference operator Cl^ with the same support as 
f{£-y ) such that 

£/Q^ = Q-J{C^). (4.7) 

We can write Qj as follows 



1=0 

where Qj '\ I — 0, ... ,k denote the coefficients in front of the shift operators E~K 
From the explicit formulas (|2.10p . (|3.4p and (13. 5p we see that for I — 1,2, ... ,k 

g^-') =0 for 7 = 0,1,..., /-I. (4.8) 

Moreover, Q^^ — ix-y^~^-k-i/X'y-i)'''-y-i and therefore equations p.6p show that 
Q^^ for 7 G Nq. This implies that the coefficients of the semi- infinite operator 
£^ are well-defined for n G No and we have 

CiQn = QnfiCn). (4.9) 

Since Pn{x) — QnPn{x), we see that equations (|2.5p and (j4.9p imply equation 
(plOall . 



5. The ^-difference equations in z 

functions of g" by 



5.1. Auxiliary facts. We define the involution on the algebra C((7") of rational 



ahcdq'^ i ^ 
From (|2.9p it is easy to see that 

In particular, In\)^n-j/2) — hence every polynomial of A„_j/2 is invariant 

under the action of In^ . Conversely, if r € <C[q"',q~^] is a Laurent polynomial of 
g", invariant under In \ then r e C[A„__,y2]- 
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Using the above, we can show that, up to a simple factor, t„ defined by equations 
p.3p - p.5|) belongs to 'C[\n-(k-i)/2\- Indeed, as we noted in Section [H we can re- 
write r„ as a fc X fc determinant whose entries are Laurent polynomials of g". 
Moreover, one can show that lit ^'^ reverses the order of the rows (i.e. it maps the 
l-ih row into the fc — / + 1 for Z = 1, 2, . . . , fc). Indeed, if 5j = a then the (Z, j) entry 
of the determinant which represents t„ is (f)^^\Xn-i+i) and we have 

/(^-i)(^W(A„_,+i)) = <P^'\K+i-k). (5.1) 
If, say, 5j = h then we need to multiply the j-th column of the Wronskian for- 
mula p.3p by 'q^'' ^n-j+i/ ^n-j+i- The (Z, j) entry of the resulting determinant is 
4'''^0(^)(A„_z+i), where 



fb pa 
ta pb 



^b, 2ipj^2ipi±l, (5.2a) 



Using the explicit formulas (I2.10p and p.4p it is easy to see that 



1 /b^'-' 



4'''' = [^-j {acq^-''+\adq'^-^+^-q)k-i {bcq^-'+\bdq^^-'+'-q)i^,. 

(5.2b) 

A straightforward computation now shows that 

Ilt-^{^^r)^^'^+^^^ (5.3) 

and the cases 5j = c or 5j = d are similar. 

Thus, Ik ' reverses the order of the rows and therefore 

This formula shows that if fc = 0, 1 mod 4, then t„ e ^[Xn-{k-i)/2]- Otherwise, 
r„ is divisible by 

K-k/2+i - K-k/2 = (1 - q)q''^^-Hq-" - abcdq"-'^) 

and the quotient belongs to C[A„_(fe_i)/2]- Therefore, we see that 

rn = ei'=^T(A„_(fe_i)/2), (5.4) 

where f is a polynomial and 

^(,)^ri iffc^O,l mod 4 ^^^^^ 

[K-k/2+i - Ki-k/2 iffc = 2,3 mod 4. 

Next we want to characterize the elements of the algebra 2l„ defined by (|3.9p and 
to show that 2l„ contains a polynomial of every degree greater than deg('f ). We can 
define the elements of 2l„ as discrete integrals of Laurent polynomials divisible by 
T„. To make this more precise, we introduce the following notation. For n,m € Z 
and for a function fs defined on Z it will be convenient to define the discrete integral 

fELm+iA ifn>m 
fsdfid{s) = S if n = m 

I - Es'ln+i if " < 
Thus 

nn—l 

fsdfid{s) = fn-\- fsdfia{s) for all m,n eZ. (5.6) 



/ 
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Note that for every polynomial r there exists r e R[i] such that 
r{K-k/2) - r{Xn-k/2~i) . 

r r = ?'(A„_(fc+i)/2j. 

Conversely, for every polynomial f there exists a polynomial r such that the above 
equation holds. Moreover, up to an additive constant, r is uniquely determined by 



r{K-k/2) ^ I (•^s-fe/2 - As_fc/2-l)^(As-(fe+l)/2)c?Md(s) 

-1 

{^s-k/2+1 - ^s-k/2)r{Xs~{k-l)/2)d^ld{s)■ 

Summarizing the above, we obtain the following characterization of 2l„. 
Proposition 5.1. For a polynomial h we have ^(A„_fc/2) £ Sin if and only if 

/"" 

h{K-k/2) = 4''^^^3(-^s-(fc+i)/2)Ts-irfA*d(s) + Co, (5.7) 



where g is a polynomial and cq is a constant. In particular, 2t„ contains a polyno- 
mial of every degree greater than deg(f). 

For the proof of (|3.10bp we need two more facts. First we formulate a discrete 
analog of a lemma due to Reach [22] , which was derived and used in [10] to establish 
the bispectral properties of specific rational solution of the Today lattice hierarchy. 
It was used more recently in [TT] |T2] within the context of orthogonal polynomials 
satisfying higher-order differential equations. 

Lemma 5.2. Let f^\fi^\-.-,fn''^^'' be functions of a discrete variable n. Fix 
rii, 71,2, ... , i^k+i G 2 o,''^d let 

^^™=E(-1)"^'^'/^'^ / /i°^Wr,(/W,..., /«,..., /r^))dMd(.), (5.8) 
j=i -^"^ 
with the usual convention that the terms with hats are omitted. Then 

/.n-1 

Wr„(/W,...,/W,i^„) 



/i«)Wr,(/ii),...,/W)dMd(s) 

-1 

xWr„(/W,...,/r^)). 



(5.9) 



Since the application in our case is very subtle and we need different elements 
of the proof, we briefly sketch it below. 



Proof of Lemma 15.^1 Note that 






Jn 

fW 
J n-1 


Jn 

f(2) 

Jn-l 


Ak+1) 
Jn 

Ak+1) 
Jn-l 




Jn-k+1 
M 
Jn-l 


f(2) 

Jn-k+1 
f(2) 
Jn-l 


Ak+1) 
Jn-k+1 
Ak+1) 
Jn-l 


Expanding the above determinant along 


k+1 

E(- 




/l^^Wr„(/(i),. 


fU) 



= 0, for every I — 0,1, . . . , k — 1. 



,k~l. (5.10) 



j=i 
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Using (|5.6p and (|5.10p we see that for / = 0, . . . , fc — 1 we have 

F„_,=^(-l)'=+i+^/^^2, / /i°)Wr,(/W,..., /«,..., /i'=+i))dA.d(.), (5.11) 
j=i -^"j 

and 

fe+i 



j=i 



(5.12) 



/rwr„(/(^\...,/r^)). 



If we phig ()5.1ip and (|5.12p in Wr„(/,1^\ . . . , f^\Fn), then most of the terms cancel 
by cohimn ehmination and we obtain (|5.9p . □ 



Remark 5.3. We hst below important corollaries from the proof of Lemma 

(i) Note that the right-hand side of (|5.9p does not depend on the integers ni,n2 ■ ■ ■ ,nk- 
Moreover, if we change Uk+i then only the value of 

" '/f'Wr,(/W,...,/W)dMd(s) 

'Mfc + l 

will change by an additive constant, which is independent of n and /A . Thus, 
instead of (15.81) we can write 



fc+l „n 

leaving the lower bounds of the integrals (sums) blank and we can fix them at 
the end appropriately. This would allow us to easily change the variable, without 
keeping track of the lower end. 

(ii) From (|5.10p it follows that for every I = —1,0,1 . . . , k — 1 we can write Fn also 
as 

Fn = / /i°^Wr,(/ii), . . . , /p), . . . , fi'+'^)d^^,is), 

and changing the variable in the discrete integral we obtain 

(iii) Suppose now that k is even. Using (ii) above with I = k/2 ~ 1 we have 

^" — Z^>> -^-l / Js+fc/2-1 "^sUs+fc/2-1' ■ • ■ ' •'s+fc/2-1' ■ • ■ ' Js+fc/2-l''"WlAJ. 

(5.13) 



Let us consider the sum consisting of the first k integrals: 

k 



^ pn 

^n ^ ^) Jn j Js+fe/2-1 "^sUs+fc/2-1' • • • ' Js+fc/2-1' • • • ' Js+fc/2-l''"f^dlij- 
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Expanding each Wronskian determinant along the last column we can wr ite F^:> 
as a sum of k terms Fn^'"^\ each one involving as integrand one of the functions 
fs+l^\ where m = — fc/2, — fc/2 + 1, . . . , k/2 — 1. We can use (IS.lOp once again, 
this time for the functions /^°'', fn \ . . . , /n*^'' (i.e. omitting fit^^''), to change s as 
follows: 



• If m > we can replace n with n — m in the upper limit of the integral, 
or equivalently, if we keep the upper limit of the integral to be n, we can 
replace s by s — m in the integrand. Thus F^'"^^ will have fs^~^^^ as 
integrand (and the integration goes up to n). 

• If m < —1 wc can replace s by s — m — 1, thus F^'"^^ will have /s*!^^'* as 
integrand (and the integration goes up to n). 



This means that we can rewrite Fn as sums of integrals, and the integrands can 
be combined in pairs (corresponding to m and (— m — 1) for to = 0, 1, ... , k/2 — 1) 

involving fs'^'^^^ and /^^i^^'. 
up to a sign, have the form 



involving fs'^'^^^ and /^^i^^'. Explicitly, we can write as a sum of terms which. 



f (0) A 1 J' f(k+l) _ f 

J s-m+k/2-l'-^s Js J: 



(0) 

s+m+fe/2 



-l"s J: 



(fc+1) 
1 



d^idis), (5.14a) 



for m = 0, 1, . . . , fc/2 — 1 and j = 1, 2, . . . , fc, where 



^(1) 

■' n~\~k/2 — 7n — l 
J n+k/2-m-2 



f(l) 



J n+k/2-m-l 
U) 

;+fc/2-m-2 



J n+k/2-m-l 
Ak) 

J n+k/2-m-2 



and 



-k/2-r. 



J n+k/2+m-l 
J n+k/2+m-2 



Jn 



A3) 

J n — k/2 — m 



(k) 



;(i) 

Jn-l 



U) 

+k/2+m-l 
A3) 

J n+k/2+m-2 



A3) 
Jn-l 



Jn 



Ak) 

-> n — k/2-~m 



Ak) 

J n+k/2+m-l 
Ak) 

J n+k/2+m-2 



Ak) 
Jn-l 



(5.14b) 



•I tl,— 



fc/24 



Aj) 

•I n,— 



fe/24 



Ak) 

•I tl, — 



fe/24 



(5.14c) 



(iv) Finally if fc is odd, we shall use (ii) with / = (fc — l)/2 and I ~ (fc — 3)/2 and 
write Fn as the average of these two sums. Then we can apply the same procedure 
that we used in (iii) to write F^^ as sums of integrals, whose integrands can be 
combined in pairs involving /i'^"''^^ and f]!^^\ 



The second important ingredient needed for the proof of (j3.10b| is the following 
results which represents an extension of Lemma 4.5 in |11) . 
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Lemma 5.4. For r{q") G C[q",q "] and £ e {a,b,c,d} there exists a q-aiuerence 
operator B^'^ such that for n G Nq we have 



[riq^^ipsix) - T{q^-^ /abcd)^i_^p,^,{x)W^{s) = B'/Cpn{x). (5.15) 



1 

Proof. Without restriction we can assume that i = a. Note first that (15.15^ for £ = a 
is equivalent to show that for every r(q") G C[g",q^"] there exists a g-difference 
operator S"^'' such that 

r{qnCPn{x) ~ r{q^-''/abcd)Cn-iPn-i{x) = Bf-[CPn{x) - C-iPn-^iix)]- (5.16) 

From the last formula it is clear that the set A of all Laurent polynomials r for 
which there exists a g-difference operator B'^-^ satisfying (|5.16p is a vector space 
over C that contains the constants. Note also that 

A„ = (g-" - 1)(1 - a6cdg"^i) G A, (5.17a) 

2 — 71 

ro(9") = 9" - ^ = -^(9""+' - abcdq^-') G A, (5.17b) 
abed abed 

by (EH) and ((2TT|) . respectively. Moreover, if r(q") G A and Riq"-) G C[g" + 
q^~'^ /abed] n A (i.e. Ri^q^) is invariant under the change g" — > q-^^^/abcd), then 
it is easy to see that rR G A, because holds with S"'''-^ = Since 

arbitrary r(g") G C[(7"',g~"'] can be written as 

r{q^')^R\q^) + R^{q")roiqn, (5-18) 
where tq is given in (|5.17bp and 

R\q^) = 1 [r((7") + riq^-'' / abed)] G C[<7" + q^-"/abcd], (5.19a) 

we see that it is enough to prove that R{q") = q" +q'^^" /abed G A. But this follows 
from equations (j5.17p and the fact that A contains the constant polynomials, since 

J?(g") - I - A I ^~^ r„(g") I + 

afocd a6c(i(l + g) 1 + 9 a6c(i(l + g) 

□ 



Remark 5.5. Note that the above proof can be used to write an explicit formula 
for the operator ^f'*" in ()5.15p as a polynomial of the operators Bz and B^ by rep- 
resenting r(g") as a sum of symmetric and skew-symmetric parts given in equation 

5.2. Construction of the g-difFerence operators in z. We are now ready to 
prove that for every h G 2t„ there exists a g-difference operator B'^ such that (j3.10bp 
holds. From Proposition 15.11 we know that, up to an additive constant, we have 

/n~l 
f'i'+l^^ 9{>'s-{k-l)/2)Tsd^J.d{s). 
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We apply Lemma [5. 21 with 

~ i-rfe+1 "+1 9[\i-{k-l)/2), 

llj=l ?n-j+l 

/l^'^ forj- = l,2,...,fc, 

If we multiply equation (|5.9I) by Xn/Iljii ?n-j+i then, using (j3.2bl) and p.Sp we 
see that the right-hand is equal to (/i(A„_fc/2) + co)pn(a;), where cq is a constant 
(independent of n and x). The goal now is to show that, if we choose appropriately 
the integers rij in Lemma 15. 2[ then there exists a difference operator Bz (with 
coefficients independent of n) such that 

F„=BzCPnix). (5.20) 

Indeed, if this equation is true, then we can pull the operator Bz in front of the 
Wronskian determinant on the left-hand side of (j5.9|) and using (j3.2bp again we will 
obtain the left-hand side of (jS.lObI) . thus completing the proof of (I3.10b[) (and, up 
to an additive constant, the operator Bz is the operator Bz)- 

Suppose first that k is even and let us write Fn as explained in Remark 15.31 
(iii). Note that the last term in the sum (15.13^ representing F„ is of the form 
(/i(A„) -I- Cg)^^p„(x), where c'q is a constant (independent of n and x), and using 
(12. 8p we can rewrite it as {h{Bz) + CQ)^^p„(a;). Thus we can consider the sum Fi''^ 
of the first k terms. It is enough to show that each term of the form (|5.14ap can 
be represented as Bz ^nPn{x) for some g-difference operator Bz- We prove this by 
using Lemma 15.41 with £ ^ Sj- Note that 

r(<z2-«/a6cd) = /W(r(g")), 

so the key point is to show that the coefficients in front of ^t'ps{x) and ^^^_iPs-iix) 
in (|5.14ap are connected by the involution Is^\ It is useful to connect the involution 
In^ to the involution In~^\ which we used earlier, by the intertwining relation 

^k/2-m-l ^ j(k-l) ^ ^g ^l) 

We have two possible cases: 5j — a ot 5j G {b, c, d}- 

Case 1: S-j = a (hence fi^^ = (A„)). Note that if we set /„ = Xn/Dj^i C-j+i^ 
then 

/(i)(/„_™+fc/2-iA^'^') - (-l)(^-i)(^-2)/2 /„+„+,/2_lA2'^ (5.22) 

where A^'^ and A^'-' are the determinants given in equations (|5.14b|) and (j5.14cl) . 
respectively. Indeed, the same arguments that we used at the beginning of Section[5] 
show that, after we multiply the j-th column of A^'-' and A^'^ by the corresponding 

factor (77^''*^^Li+i/Cn-i+i) from / for every i such that Si ^ a, then In^ will reverse 
the order of the rows in these two determinants (and formula (|5.2ip shows that the 
arguments work in the same way with /i*^ there replaced by In"^ here.). Since 
Ali^ are (fc — 1) x (fc — 1) determinants, we obtain (I5.22[) . Note also that 

-^n^'(5(An-m-l/2)) = ff(A„+m-l/2) 

and 

Kl). (fe+2) ^ (,^(k-l){k-2)/2 (k+2) 
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The last equation, follows easily from (|5.5|) by considering separately the two pos- 
sible cases k = mod 4 and k = 2 mod 4. Since 

/i°) = /„6i';+^)5(A„-(fc-i)/2), 

we see that Lemma with £ — a can be applied for the integral in (|5.14a[) and we 
can rewrite (I5.14ap as 

for some g-difFerence operator B'^. We can now commute B'^ and </>'■-'•' (A„) and use 
(|2.8p to rewrite the above equation as 

B'J^^HB.)CPnix), 
completing the proof in this case. 

Case 2: Sj G {b,c,d}. Let us fix 5j ~ b (the cases Sj — c or 5j — d are simi- 
lar). Then fn"^ — i,n4>^^\^n) I S,n- Notc that the j-th columns of the determinants 
A^^^ and A^'-' are removed and therefore the main difference between this and the 
previous case is that in fn^ we have the extra factor '7^'-'^^_j+i/'?^_j+i- Using the 
arguments in case 1, we see that the integral in equation (I5.14ap can be rewritten 
as 

's—m+k/2 — l^s—m+k/2—i / ^\^n / \ 



fra 



ls+m+k/2-l^s+ni+k/2-j ^(i). 



fa 

^s+m+k/2-j 



s—m 



where in the last equality we used (|5.2aP . In order to apply Lemma 15.41 with i — b 
we need to check that 

T{l}(b;k/2-m,j n _ b;k/2+m+l,j 
\'^n-m+k/2-l) '^n+m+k/2-1^ 

which follows from equations (|5.3p and (|5.2ip . Thus we can construct a g-difference 
operator B'^ such that the expression in (I5.14al) can be rewritten as 

fi'^KC.Pn{x) = B'A'^^\Xn)^,Pn{x) = B'^^^^^B CnPn{x) , 

completing the proof when k is even. The case when k is odd follows along the 
same lines, using the representation of Fn described in Remark l5.3l (iv*). □ 

Remark 5.6. Using Remark 15. 5[ we can write an explicit formula for the q- 
difference operator B'^ in equation (I3.10b[) in terms of the operators Bz, B^, £ € 
{a, b, c, d}. This might be useful for specific examples or small values of fc, but the 
formula becomes very involved in general. 

6. BiSPECTRAL EXTENSIONS WHICH ARE ORTHOGONAL WITH RESPECT TO A 

MEASURE ON THE REAL LINE 

In this section we assume (for simplicity) that the parameters a, 6, c, d are such 
that max(|a|, |c|, |d|) < 1. Then the Askey- Wilson polynomials {p„(a;; a, 6, c, d)}5^o 
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with 



are mutually orthogonal with respect to the measure ^a,b,c,d on [—1,1], defined by 

w{x;a,b,c,d) 

diia,b.c.d = — — , „ ax, (6.1a) 

ZTTV 1 — X'' 

w{x; a, b, c, d) = , , , ,' ' , °° , (6.1b) 

(az, a/z, oz, o/z, cz, c/z, dz, d/z; q)oc 

see [a Theorem 2.2, p. 11]. 

Next we explain how the techniques developed in the paper can be used to obtain 
extensions of the Askey- Wilson polynomials which eigenfunctions of higher-order 
g-difference operators and are mutually orthogonal with respect to a measure on R. 
This provides a new derivation of the constructions in our joint paper with L. Haine 
[5] . Applying Theorem 13.11 to this situation, we obtain g-difference operators of 
lower-order than the ones constructed from the general theory in [9 . 

The key point is to choose the polynomials (A-y) so that the space W de- 
fined by p.7|) is C~f invariant (and therefore, the polynomials Pn{x) will satisfy 
a three-term recurrence relation). The simplest possibility is to have 
eigenfunctions for the operator which in view of equation (|4.ip is equivalent to 



pick the polynomial ^^■'■'(A-y) to be an eigenfunction for the operator ct^' for every 



j. Note that Proposition 14.21 alreadv gives us polynomial eigenfunctions for the 
operators Ci^ where I S {a, 6, c, d} and therefore one natural choice is, for every 
J = 1, 2, . . . , fc, to take 

0(^)(A,)=5j(A,), (6.2) 

where kj E Nq. However, the corresponding extensions of the Askey- Wilson poly- 
nomials will depend only on the four parameters (a, b, c, d) and will not contain as 
a limit the Krall-Jacobi or even the Krall-Laguerre polynomials. 

In order to get more interesting examples, we fix ^ e {a, 6, c, d} and we try to 
find special values of the parameters a, b, c, d for which the difference operator £^ 
has two dimensional polynomial eigenspaces (for specific eigenvalues). This would 
allow us to pick (l)^^\X-y) in the corresponding eigenspaces (and they will depend 
on new free parameters). 

Recall that the Askey- Wilson polynomials Pnix) — Pnix] a, b, c, d) are symmetric 
in the parameters a, 6, c, d and therefore we can rewrite formula (|2.2p by exchanging 
the roles of a and c?, i.e. we have 



{ad,bd,cA;q)n 

Pni^) = J. 4<i 



q ^,abcdq^ ^,dz,dz ^ 

ad, bd, cd ' ^' ^ 



(6.3) 



The series above clearly terminates for all 71 G C if we fix z = = dg™ where 
m e Nq. Set Xm — ■^{zm + J-)- Equations (|2.5p and (|4.ip show that 



9 {abcd/q,ad;q)n . 

UPn{Xm) = 7-7^ 77 T 403 

(ad)" {bc,q;q)n 



q ^,abcdq" ^,dz,n,q ™ 
ad, bd, cd ' 



(6.4) 



is an eigenfunction for the operator £° with eigenvalue 2xrn = Zm + j-- Moreover, 
since the 4^3 series terminates, it is easy to see that the meromorphic extension of 
the expression in (j6.4p is an eigenfunction for the bi-infinite extension (because 
after we cancel the exponents and the g-shifted factorials in front of the 4^3 series, 
the eigenfunction equation can be rewritten as a polynomial equation in q"^ , which 
if true for j — n e Nq, must be true for all 7 G C). Note also that, if a = dq" 
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where a G N, the eigenvalue Zm + j- — dq"^ + above, and + = 
dq^3+a-i _|_ in. (|4.5p coincide when m = —j + a— 1. Moreover, if ad = (7"+' 

where / G N then the exponents and the q-shifted factorials multiplying the 4^3 
series in (j6.4l) can be rewritten as a Laurent polynomial in as follows 

g" {abcd/q,ad;q)n 1 (focg", 9"+^; 



(ad)" (&c,Q;g)„ (foe, g; g)„+i_i (g»)a+i-i 

It is not hard to see that the above expression is a polynomial of A„ (since it is 
invariant under the action of the involution In"^) and that the eigenfunctions in 
(|4.4p and (|6.4p are independent. 

Summarizing the above comments, we obtain a different proof of Proposition 5.2 
ini. 

Proposition 6.1. Assume that d — eq''^^ and a = dq°' = £(7"+'/^, where e — ±1, 
and a, / G N. For to = 0, 1, . . . , a — 1, the polynomials 



U^'^iX^) = 403 



g-"+"+Va2, g-^+i, l/abcdq'<-^ . 



q^/ab.,q^/ac,q^/ad ' ^' * 



2,m/^ N (focgT,gT+^g)a+;-i 
^ 1-^7) - /_.A^j-;-i 403 



ad, M, cd ' ^ 



form a fundamental set of solutions for the eigenvalue problem 



(6.5a) 
(6.5b) 



C^u^ = 2xmU^, where x,„ = -(g™+'/^ + 9"™"'/^). (6.5c) 

The polynomials u^'™, u^'™ above correspond to s.y(z„i) and r^{zm) in [9J page 302]. 
The connection between u^'™- and s^(z„i) can be seen by applying Sears' transfor- 
mation formula page 49, formula (2.10.4)]. 

If the parameters a and d satisfy the conditions in Proposition 16. 1[ then we can 
pick 5j = a and 0(-'^(A^) G spanju^^^^ (A.y), -u^'™^ (•^7)}i for some ruj G {0, 1, ... , a— 
1}. Note that span{0(-') (A^)} depends on one free parameter. If we make a similar 
choice for every j — l,2,...,fc, then clearly C-^W C W (since each ■^^''^ is an 
eigenfunction for C-y). We shall assume that the polynomials ^^•'^(A-y) are generic 
so that the conditions p.6p hold. Therefore 21^ = C[z + I/2] and the polynomials 
{Pn(a;)}J^o ^ill satisfy the three-term recurrence relation 

CnPn{x) = 2xpn{x), whcrc £„ := 

It is well knowrQ that the intertwining relation (14. 7p means that the operator 
can be obtained from C-y by a sequence of fc elementary Darboux transformations 
at the points ■ For the semi-infinite operators, this gives 

^71 — ^milci + ^ Qn ^ ^ii ^ ^ ^milci -|- ^ Pfi ~ ^7712^^ ^ Qn ^ ^ ' ' ' 

Ci'-'^ = a;,„,_,Id + Qi^-'Ti'-'^ = x„7jd + )Ql'^-) (6.6) 

where Pn '^ ~ Pn'' En + p"'-'Id are first-order forward difference operators, QiP = 
q"^^Id + q^^'^E~^ are first-order backward difference operators, and 



^See O Proposition 6.1], although this statement can be found aheady in |23l pp. 14-19] ir 
the case of differential operators. 
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If we define polynomials Pn\x) — Qn^Pn where pn\x) — p„(x) are the 

Askey- Wilson polynomials, then pn{x) — pit\x) are the polynomials constructed 
in Section 131 Using (|6.6p . it is not hard to prove inductiveljQ that the polynomials 
{pn\x)}'^^Pi are orthogonal with respect to a measure of the form 

7 (i) d^a,b,c,d V"^ c/ \ 

where iJa,b,c,d is the Askey- Wilson measure in (j6.ip and the constants 1^1,..., Vj 
are in one-to-one correspondence with the free parameters that specify the one 
dimensional subspace spanlt/)^-'^ (A^)} of spanju^'^^ (A^), u^'™^ (A^)} at each step. 
Putting all this together, we obtain the following theorem (c.f. Theorem 6.2 in [9 ). 

Theorem 6.2. Suppose that d = eq^^^ and a — £g"+'/^ , where e — ±1, and a,l £N 
and let k < a. 

(I) For each j = 1, 2, . . . , fc pick Sj — a and 

(t)'^^H\-y) € span{Mi'™^(A^),w^''"'(A^)} (6.7) 

with rUj a {0, 1, . . . , a — 1} and rrij ^ mi for j ^ i. Then for generic {(f>^^^ (''^7)} as 
in (|6.7p . the polynomials Pn{x) defined by (|3.ip - (l3.2al) are orthogonal with respect 
to a complex measure of the form 

df^ = ^f^P^^+il'^A^-^^.) ^'ih = £(g™.+V2+^-™,-;/2)^ (g 8) 

(II) Conversely, ifnij G {0, 1, . . . ,a — 1} are distinct numbers, then for generic com- 
plex numbers vi, . . . ,Vk we can use p.ip - (l3.2al) to construct polynomials {Pn{x)}^^Q 
orthogonal with respect to the complex measure in (j6.8p by picking Sj = a and an 
appropriate ip'^^^Xy) as in (j6.7p for every j — 1,2, . . . ,k. 

We can use Theorem 13.11 to deduce that the polynomials Pn{x), which are or- 
thogonal with respect to the measure in (|6.8p , are eigenfunctions of the g-difFerence 
operators in the algebra D^. In the above theorem, we can replace a and d with 
any two of the parameters (a, b, c, d). In particular, if two of the parameters satisfy 
the above conditions with e = 1, and the other two satisfy analogous condition with 
e = — 1, we can add point masses on both sides outside of the interval [—1,1]. More 
precisely, suppose that b — eg^+*/^, c — eq*/^, where G N and e = ±1. Then 
for m = 0, 1, . . . , /3 — 1 we can define functions v^'™{\^) and v'^''^{X'y) by formulas 
(j6.5ap and ()6.5bp . respectively, by replacing (a, b, c, d) with (6, a, d,c), a with f3, and 
I with t. This leads to the following theorem (cf. Corollary 6.6 in [5]). 

Theorem 6.3. Suppose that a = qa+i/^^ ^ _ _^,3+t/2^ ^ ^ —q*'^^, d = q'/^, where 
a, I3,l,t Cz N. Let ki, k2 (z N be such that fci < a, fc2 < /3 and set k — ki + k2. 
(I) For each j = 1, 2, . . . , fci pick Sj = a and 

0(^'HA^) e span{ui^"^'(A^),w2.'",'(A^)} (6.9) 

with mj G {0, 1, . . . , a — 1} and mj 7^ mj for j 7^ i, and for each j — ki + 1, . . . ,k 
pick Sj = b and 

<j)^^HX^) G span{t;i'™?(A^),w2^"?(A^)} (6.10) 



^See Theorem 2 in [7]. 
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with iTij e {0, 1, . . . , /3 — 1} and mj ^ nif for j ^ i. Then for generic (A-y)} as 
in ()6.9p - ()6.10p . the polynomials Pn{x) defined by (|3.ip - p.2ap are orthogonal with 
respect to a complex measure of the form 

dfi = — ^1 ^ ' ^-9'^+ ' h v]5(x — x^i) + v'^Six — x^a), with 



.11) 



(II) Conversely, suppose that nij G {0, 1, . • . , a— 1}; ?7i| £ {0, 1, • . . , [3 — 1} and 7^ 
m[ /or j 7^ i, r = 1, 2. T/ien /or generic complex numbers v\, . . . ^v], , i/^, . . . , 
we can use (I3.ip - (l3.2ap to construct polynomials {Pn(,x)\^^Q orthogonal with respect 
to the complex measure in (j6.11l) by picking 

• Sj — a and an appropriate (A-y) as in (16. 9p /or every j — 1,2, ... , ki, 

• Sj ~ b and an appropriate (f)^^"^ (A^) as in (I6.10p /or every j — ki + l,...,k. 

Remark 6.4. As we noted earlier, we can use ()6.2p to add Dirac measures without 
free parameters. We can mix this with the constructions in Theorems 16.21 and 
16.31 to get even more comphcated examples of orthogonal polynomials satisfying 
higher-order g-difFerence equations. 

Remark 6.5. Recall that the substitution 

a = q"+i/2, 6 = -/+i/2, c--gi/2^ d^q'/\ 

in the Askey- Wilson polynomials leads to q-extensions of the Jacobi polynomials, 
see [211 . In particular, this explains how the polynomials {p„(a;)}^Q in Theorems 
16.21 and 1 6.31 can be considered as extensions of the Krall- Jacobi polynomials. 



Example 6.6. Consider the simplest application of Theorem 16.21 by taking e = 
a = / = 1, or equivalently a = q'^^'^ , d = q^^^ , b, c arbitrary free parameters. Then 
k — 1, mi — and 

0(l)(A^) = u^'^iX^) + D,u^'°{X^) = 1 + - bcq')i_^ - ^ 

where is a free parameter. Then r„ — (/''■^■'(A„) is a linear polynomial in A„ and 
therefore the algebra 2l„ will be generated by two polynomials of degrees 2 and 3, 
respectively. This means that the corresponding algebra 'Dz will be generated by 
two g-difference operators of orders 4 and 6, respectively. All this follows immedi- 
ately from Theorem 13.11 while the operator of minimal order that we can obtain 
from the general theory in ^ is of order 12, see [51 Section 7]. 

This example can be considered as a g-extension of the Krall- Jacobi polynomials 
[18) . The polynomials p„ (x) are orthogonal with respect to a measure of the form 

dp,„3/2 h,c.q'/2 , J., \ 1/ 1/2 , -l/2\ 

d^ = — 7 — ^^-T h I'lSix - xq) with xq = -(g ' + q ' )■ 

(x - ) 2 

Note that x — xq — ~ 2gV^ ^ ^ zq^l'^)(l — q^l"^ jz), and therefore fi can also be 
rewritten as 

dp. = -'2^dlJ,qi/2^i,^^^gl/2 + ViS{x - Xq). 

The connection between v'l and vi can be derived from the orthogonality condition 
Piix) ± 1. 
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Remark 6.7. Note that the discrete parts of the measures in Theorenis l6.2l and l6.3l 
are supported at specific points of the form ±^{q^^^ + q^^^^), s e N, outside of the 
interval [—1, 1]. They accumulate to the boundary ±1 only in the limit g — > 1. In 
the recent paper [1,, the authors consider different extensions of the Askey- Wilson 
polynomials, not related to the bispectral problem, by adding mass points at the 
end points ±1 of the interval [—1, 1]. 
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